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CHARACTERIZATION  OF  IONOSPHERE  WA  VEGUIDE  PROPAGATION 
BY  MONITORING  HAARP  HE  TRANSMISSIONS  IN  ANTARCTICA 

Project  manager:  Yuri  M.  Yampolski 
1.  Introduction 

The  nonmonotonic  run  of  the  height  profile  of  the  ionospheric  electron  density  associated 
with  the  multilayered  structure  of  the  ionosphere  creates  preconditions  for  formation  of  interlayer 
duct  channels.  These  natural  structures  are  capable  of  supporting  the  radio  wave  energy  transfer  to 
long  distances  (thousands  and  tens  of  thousands  of  kilometers)  with  low  losses.  The  valley 
between  the  E-  and  F-layer  maxima  seems  to  be  energetically  optimal  from  this  point  of  view.  The 
characteristic  boundaries  of  such  a  duct  vary  between  120  and  150  km  on  the  bottom  and  from 
200  to  250  km  on  the  top,  depending  on  the  wave  frequency  and  geophysical  conditions 
(geographic  coordinates,  sunlit  conditions,  season,  geomagnetic  activity  etc.).  At  HF  waves  the 
duct  represents  a  multimode  waveguiding  structure  since  the  wavelength  X  is  much  shorter  than 
the  waveguide  size  h.  The  ionospheric  duct  is  characterized  by  a  small  amount  of  losses  owing  to 
the  comparatively  low  collision  frequency  of  electrons  Ve  whose  value  is  two  or  three  orders  of 
magnitude  lower  than  at  heights  of  the  main  absorbing  D-region  of  the  ionosphere.  This  is  a 
reason  why  the  attenuation  rate  of  HF  signals  in  the  interlayer  duct  is  dozen  of  times  smaller  than 
for  ordinary  multi-hop  propagation  paths. 

The  major  problem  of  using  the  ionospheric  duct  in  practice  is  the  difficulty  of  exciting  the 
channel  from  the  Earth’s  surface  and  extracting  energy  for  a  ground-based  consumer.  The  reason 
is  that  the  E-layer  plays  the  role  of  a  barrier  which  shades  the  above  layers  from  the  HF  emission 
radiated  by  a  ground-based  transmitter  at  frequencies  below  the  critical  frequency  of  the  E-region. 
The  signal  input  in  and  extraction  from  the  interlayer  duct  channel  is  possible  only  in  the  presence 
of  considerable  horizontal  gradients  of  the  electron  density  or  intense  scattering  irregularities.  For 
the  known  geometry  of  the  propagation  path  it  is  possible  to  select  the  most  favorable  conditions 
for  feeding-in  (feeding-out)  the  waveguide  proceeding  from  the  variable  sunlit  conditions  and 
hence,  the  predictable  horizontal  gradients  of  the  electron  density  distribution.  Such  conditions 
can  be  mainly  realized  for  rather  long-range  radio  paths  of  meridional  orientation  close  to  the 
sunrise  or  sunset  at  ionospheric  heights.  Pointing  the  transmitting  and  receiving  antennas  in 
certain  oblique  directions  and  changing  the  radiation  frequency  it  is  possible  to  arrange  a  scenario 
where  the  bottom  wall  of  the  waveguide  (E-region)  would  exist  along  the  entire  radio  path,  except 
the  ascending  (descending)  section  of  the  trajectory  in  a  given  vicinity  of  the  transmitting 
(receiving)  site.  Such  situations  are  well  known  to  occur  in  practice  of  radio  broadcast  and 
communication  when  the  solar  terminator  passes  almost  simultaneously  the  transmitter  and 
receiver  sites.  This  can  be  accompanied  by  a  considerable  increase  in  the  received  signal  level. 

The  effects  of  the  energy  input  and  output  owing  to  radio  wave  scattering  by  intense 
natural  ionospheric  irregularities  are  of  unpredictable  character.  Their  appearance  is  determined  by 
sporadic  factors  and  depends,  specifically,  on  the  enhanced  turbulization  of  the  ionopsheric 
plasma  produced  by  various  kinds  of  geophysical  perturbations.  Especially  frequent  such  effect 
are  at  high  latitudes  which  are  to  a  greater  extent  subjected  to  the  influence  of  magnetic  and 
corpuscular  perturbations  stimulated  by  solar  wind  flux  variations.  It  is  evident  that  predictable 
input  (output)  of  the  useful  signal  into  (from)  the  duct  using  the  effects  of  scattering  by  natural 
structures  of  the  kind  is  impossible.  In  addition,  the  life  time  of  the  necessary  irregularities  is 
unpredictable  as  well  and  varied,  as  a  rule,  from  a  few  seconds  to  dozen  of  minutes. 
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Extra  chances  for  exciting  the  interlayer  duct  can  be  provided  by  active  experiments  on  ionosphere 
modification  which  allow  controlling  the  space  and  time  structure  of  the  ionospheric  turbulence. 
Primarily  this  concerns  the  ionospheric  heating  by  the  powerful  HE  emission  from  a  ground-based 
transmitter.  As  is  known,  the  major  result  of  the  nonlinear  interaction  between  the  incident  wave 
and  ionospheric  plasma  is  generation  of  the  artificial  ionospheric  turbulence  (AIT),  whose  space- 
and-time  spectrum  can  be  controlled  using  special  heating  regimes.  As  a  rule,  the  horizontal  scale 
of  the  AlT-containing  region  is  determined  by  the  antenna  pattern  width  of  the  heater  and  can  be 
estimated  to  lie  between  100  and  200  km,  depending  on  the  modified  region  altitude  and  the 
heating  frequency.  It  is  evident  that  the  stimulated  irregularities  are  of  a  stochastic  character, 
however  a  number  of  their  morphological  features  are  a  priori  known.  These  are  the  location; 
times  of  appearance,  existence  and  relaxation;  spatial  anisotropy  due  to  the  magnetized  plasma; 
external  scale-size  of  the  turbulence;  and  even  spatial  spectrum  in  the  case  of  special  heating 
modes.  With  allowance  for  such  deterministic  features,  the  AIT  can  be  used  for  exciting  the 
interlayer  channel.  The  output  of  the  signal  from  the  waveguide  can  be  provided  by  horizontal 
gradients  produced  by  the  solar  terminator  near  the  receive  site  of  a  super  long-range  path. 

A  pioneering  feature  of  the  developed  theoretical  model  is  the  account  of  the  regular 
ionospheric  refraction.  The  aspect-sensitive  contours  in  the  ionosphere  and  on  the  Earth’s  surface 
have  been  calculated  for  all  the  transmitting  and  receiving  sites  for  the  current  ionospheric 
conditions.  During  the  measurement  campaigns  with  the  use  of  radiations  from  the  powerful 
heaters  two  other  effects  have  been  occasionally  revealed,  which  are  the  combination  (Brillouin) 
signal  scattering  by  the  artificially  stimulated  plasma  turbulence  and  excitation  of  the  second 
harmonic  of  the  powerful  radiation  in  the  HE -modified  ionosphere.  As  planned,  a  software 
package  has  been  developed  for  the  remote  control  of  the  receiving  complexes  and  visual 
representation  of  the  measurement  results  in  real-time  through  the  internet. 

Taking  into  account  that  the  most  powerful  heaters  HAARP  and  EISCAT  are  located  close 
to  the  northern  polar  regions  (Alaska  and  Northern  Scandinavia,  respectively),  the  most  promising 
position  to  observe  the  waveguide  propagation  effects  of  HE  signals  to  super-long  distances  seems 
to  be  Antarctic.  An  additional  reason  in  favor  of  this  region  is  the  “purity”  of  the  electromagnetic 
climate  because  of  the  absence  of  intense  man-caused  interference  and  local  thunderstorm  activity. 
The  first  successful  experiments  on  registration  of  signals  from  a  powerful  heating  facility  in 
Antarctic  region  were  performed  by  the  Institute  of  Radio  Astronomy  in  2002.  The  results 
suggested  that  the  waveguide  propagation  can  apparently  be  realized  Zalizovski  et  al,  [2009].  The 
probe  signal  was  the  proper  emission  from  the  EISCAT  heater  which  was  scattered  by  AIT  inside 
the  interlayer  duct.  The  receiving  site  was  deployed  in  Antarctica  at  the  Ukrainian  base 
“Akademik  Vernadsky”.  In  a  number  of  cases,  along  with  the  usual  quite  stable  component 
associated  with  multi-hop  propagation  through  the  sidelobes  of  the  heater  antenna,  rather  strong 
signals  were  reliably  detected  whose  spectra  allowed  suggesting  that  they  were  scattered  by  AIT 
produced  by  the  heater  emission  itself  Most  frequently  the  scattered  component  was  observed 
close  to  the  time  of  solar  terminator  passage  through  the  radio  path.  The  effect  has  been  called  the 
“self-scattering”  of  powerful  HE  signals  by  the  ionospheric  turbulence  stimulated  by  the  same 
signals  Zalizovski  et  al,  [2009].  However  these  were  trial  episodic  experiments,  which  have  not 
allowed  investigating  in  detail  various  mechanisms  of  radio  signal  propagation. 

The  present  Project  was  aimed  at  comprehensive  studying  the  super  long-range 
propagation  mechanisms  through  monitoring  the  HAARP  and  EISCAT  heaters  emissions  with  the 
use  of  a  specially  constructed  automated  internet-controllable  HE  receiver.  Possible  receiving  site 
should  been  deployed  in  Antarctica  at  the  Ukrainian  Antarctic  station  (UAS)  “Akademik 
Vernadsky”  (65®  15’  S,  64®  16’  W).  The  distances  between  HARRP  and  EISCAT  heaters  and 
Antarctic  position  are  15  600  km  and  16  100  km  correspondently. 
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2.  Mobile  HF  receiving  system 

For  the  Stage  1  of  the  Projeet  the  mobile  HF  receiving  complex  was  created  and  installed  at 
the  Low  Frequency  Observatory  of  IRA  NASU(LFO)  for  tests  and  continuous  measurements.  The 
receiving  facility  has  been  budett  around  the  WR-G313i  digital  receiver  of  the  WiNRADiO 
communications  (Australia,  www.winradio.com.au).  which  has  shown  good  performance.  The 
receiver  should  be  designed  as  a  PCI  card  to  be  inserted  into  and  extender  slot  of  the  PC 
mainboard.  A  general  view  of  the  receiver  is  shown  in  Fig.  2.1. 


Fig.  2.1.  General  view  of  the  WR-G3 13i  receiver 
The  basic  performance  specifications  of  the  receiver  are  listed  in  Table  2.1. 


Table  2.1  Performance  specifications  of  the  WR-G313i  receiver 


Receiver  type 

Superheterodyne  receiver  with  a  signal  processor 

Frequency  range 

9  kHz  to-  30  MHz 

Frequency  adjustment  accuracy 

1  Hz 

Dynamic  range 

95  dB  (1 12  dB  with  an  input  attenuator) 

Sensitivity 

0.05  pV 

Frequency  band 

Variable  between  1  and  1500  Hz  with  1  Hz  step 

Intermediate  frequencies 

IFl:  45  MHz 

IF2:  16  kHz  (adjustable  between  12  and  22  kHz) 

IFl  bandpass  filter 

15  kHz  crystal  filter 

Antenna  input 

50  Q 

Form  factor 

2/3  PCI  card,  PCI  2.2  compatible 

Weight 

330  g 
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o 

To  provide  for  the  required  frequency  instability  of  the  receiver  AF/F  no  worse  than  10'  ,  a 
temperature-compensated  compact-size  crystal  oscillator  developed  at  the  IRA  NASU  will  be 
used.  The  oscillator  is  designed  as  a  standard  5”  PC  device  to  be  inserted  in  a  computer  (Fig.  2.2). 


Fig.  2.2.  Compact-size  temperature-compensated  crystal  oscillator 

An  active  loop  antenna  developed  at  the  IRA  NASU  is  to  be  used  as  to  receive  the  signals.  The 
antenna  provides  for  a  high  sensitivity  within  the  frequency  range  3  to  15  MHz.  The  design  of  the 
loop  antenna  with  the  pre-amplifier  is  shown  in  Fig.  2.3 


Fig.  2.3.  Active  loop  antenna  (1,)  broadband  pre-amplifier  (2),  and  fixing  parts  (3) 

A  special  software  package  developed  at  the  IRA  NASU  allows  remotely  controlling  the  receiver 
and  real-time  data  reading  through  the  Internet.  To  detect  the  useful  signal,  Doppler  spectra  of  the 
received  emission  are  computed  within  a  narrow  frequency  range  (about  100  Hz)  with  a  spectral 
resolution  ~0.1  Hz.  Fig.  2.4  shows  the  user  window  interface  intended  for  frequency  tuning  of  the 
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receiver  through  the  Internet  and  visual  representation  of  the  current  measurement  data  in  the  form 
of  the  signal  envelope  (a)  and  instantaneous  (b)  and  dynamic  (c)  spectra. 


Fig.  2.4.  User  window  interface  of  the  receiver  control  and  data  reading 


A  general  view  of  the  receiver  with  a  control  computer  is  shown  in  Fig.  2.5. 


Fig.  2.5.  General  view  of  the  receiving  complex:  1  -  WR-G313i  receiver,  2  -  reference  oscillator, 

and  3  -  antenna  cables 

The  similar  receiving  complex  was  created  for  installation  in  Antarctica  on  the  Stage  2  of  the 
Project.  On  March  2013  this  system  deployed  at  the  Ukrainian  Antarctic  station  “Akademik 
Vernadsky”  and  operates  continuously  registering  now  emissions  from  HF  broadcasting  stations 
and  signals  from  the  HAARP  and  EISCAT  heaters. 


Distribution  A:  Approved  for  pubiic  releas;  distribution  is  unlimited. 


4/5/2015 


STCU 


PARTNER  PROJECT  P-524  EINAL  REPORT 


PAGE  6 


3.  HF  wave  scattering  by  field-aligned  plasma  irregularities  considering  refraction  in  the 

ionosphere 

The  effect  of  ionospheric  refraction  on  the  scattering  of  high  frequency  (HF)  signals  by  random 
field-aligned  irregularities  in  the  upper  ionosphere  was  analyzed  in  a  frame  of  a  theoretical  part  of 
the  Project.  The  possibility  of  excitation  of  the  ionospheric  interlayer  waveguide  by  the  aspect- 
sensitive  scattered  HF  signals  is  analyzed  in  detail  for  the  specific  conditions  of  the  HF  heating 
experiments.  Aspect-sensitive  scattering  of  high-frequency  (HF)  (3-30  MHz)  and  very  high- 
frequency  (VHF)  (30-300  MHz)  waves  by  field-aligned  irregularities  of  the  ionospheric  plasma 
are  widely  used  for  diagnostics  of  turbulent  processes  in  the  near-Earth  plasma  of  both  natural 
[Lyon,  1965;  Bourdillon  et  ah,  1995;  Bezrodny  et  ah,  1997]  and  artificial  [Djuth  et  ah,  2006; 
Hysell,  2008;  Yampolski  et  ah,  1997;  Koloskov  et  ah,  2002]  origins.  Theoretical  premises  for  this 
method  of  diagnosing  the  inhomogeneous  structure  of  the  ionosphere  were  developed  within  the 
single-scattering  approximation  [Rytov  et  ah,  1987].  As  a  rule,  in  the  relevant  studies  the 
condition  »  co]^  is  assumed,  where  co  is  the  cyclic  frequency  of  the  sounding  signal, 

is  the  critical  frequency  of  the  ionospheric  layer  (i.e.,  the  maximum  plasma  frequency  ), 

and  Wjj  the  electron  gyrofrequency.  This  assumption  makes  it  possible  to  ignore  the  effects  of  the 
regular  ionospheric  refraction  and  background  geomagnetic  field.  Consequently,  the  wave 
trajectories  are  approximated  by  straight  lines,  and  ionospheric  plasma  is  assumed  isotropic. 

However,  to  be  precise,  in  order  to  ignore  the  effects  of  refraction  stronger  conditions  have 
to  be  satisfied,  which  can  be  derived  as  follows.  The  equation  for  the  incident  ray  trajectory  in  a 
plane-stratified  ionosphere  can  be  written  as  [Kravtsov  and  Orlov,  1990]: 

x(z)  =  sin6() 

where  is  the  angle  of  incidence  at  the  lower  boundary  of  the  ionospheric  layer,  Zq  ,  counted 

from  the  vertical.  Obviously,  in  order  for  a  wave  to  penetrate  through  the  ionosphere  (with  no 
reflection  point  at  which  the  ray  curvature  would  be  the  greatest),  the  following  condition  must  be 
satisfied:  a>  >  /cos6>q  ;  while  in  order  to  approximate  the  ray  trajectories  by  straight  lines,  an 

even  stronger  condition  should  hold:  co  »  /  cos^q  .  Here  is  an  example  from  practice.  During 

most  of  the  experiments  on  HF  signal  scattering  from  artificial  ionospheric  turbulence  (AIT) 
produced  by  the  Sura  HF  heater  (Nizhny  Novgorod,  Russia)  [Yampolski  et  ah,  1997;  Myasnikov 
et  ah,  2001;  Koloskov  et  ah,  2002],  the  geometry  was  such  that  »  70°  (for  the  probe  signal)  and 

fcr  -  ^cr  /(2:t)  >  4.5  MHz.  Therefore,  in  order  to  be  able  to  ignore  the  effects  of  refraction,  the 

probe  signal  frequency  should  be  much  higher  than  ~13  MHz,  which  was  rarely  the  case. 
Therefore,  for  correct  analysis  and  interpretation  of  experiments  of  the  kind,  the  refraction  effects 
need  be  taken  into  account. 

It  is  especially  important  to  account  for  the  ionospheric  refraction  in  analyzing  the 
scattering  of  powerful  HF  emissions  on  the  ionospheric  irregularities  produced  by  the  same 
radiation.  The  fact  that  HF  pump  wave  scatters  on  ionospheric  irregularities  is  rather  trivial  and 
has  been  observed  a  long  time  ago  (Belikovich,  et  ah, 1975;  Erukhimov  et  ah,  1980).  It  was  not 
clear,  however,  whether  the  scatter  signals  can  travel  over  the  large  distances  away  from  the 
heater,  for  example,  through  the  ionospheric  waveguide.  Such  effect  was  apparently  first 
experimentally  discovered  by  Zalizovsky  et  al.  [2009]  using  the  EISCAT  HF  heater  (Tromso, 
Norway)  and  was  given  the  name  “self-scattering  effecf’.  In  their  experiment  the  EISCAT 
transmission  was  monitored  at  three  greatly  dispersed  receiving  sites:  at  the  Ukrainian  Antarctic 
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Station  “Akademik  Vernadsky”  (UAS);  at  the  Radio  Astronomical  Observatory  of  the  IRANASU 
(RAO)  near  Kharkov,  Ukraine;  and  near  Sankt-Petersburg,  Russia  (SPB).  Typically,  the  received 
signals  contained  two  characteristic  spectral  components.  One  component  was  a  relatively  narrow 
line  with  insignificant  Doppler  shift  variations  characteristic  of  the  sky-wave  HF  propagation  in 
middle  latitudes.  The  spectrum  of  the  other  component  was  much  broader  reminiscent  of  that  of  a 
signal  scattered  at  frequencies  above  the  maximum  usable  frequency  (MUF).  The  broadband 
(scattered)  spectral  component  exhibited  strong  frequency  fluctuations  with  amplitude  rate 
occasionally  exceeding  10  Hz.  However,  the  most  notable  observational  fact  was  that  the  Doppler 
frequency  shifts  and  intensities  of  this  component  varied  practically  synchronously  at  all  three 
receive  sites  with  nearly  the  same  oscillation  rate.  The  authors  concluded  that  the  observed 
Doppler  shifts  of  the  broadband  spectral  component  were  likely  induced  within  a  wave  trajectory 
segment  that  was  common  to  all  the  propagation  paths,  i.e.,  where  the  HF  pump  wave  travels  from 
the  heater  to  the  scattering  region  with  the  AIT.  Simultaneous  observations  by  the  EISCAT 
incoherent  scatter  radar  showing  strong  variations  in  the  electron  density  above  the  HF  heater  are 
supporting  the  suggested  mechanism.  However,  the  exact  mechanism  of  the  super  long  range 
propagation  of  the  HF  heater  signals  (the  EISCAT-UAS  propagation  distance  was  ~1 6,300  km) 
still  has  not  been  determined,  and  it  is  quite  possible  that  ionospheric  refraction  may  play  an 
important  role  since  the  experiments  at  the  EISCAT  location  were  conducted  at  frequencies  below 
the  critical  frequency  of  the  ionosphere  {co<co^ ).  Galushko  et  al.  [2008]  had  also  observed  the 

self-scattering  effect  by  monitoring  the  HAARP  transmission  at  several  remote  sites  in  the  USA, 
Europe,  and  Arctic  using  digital  Doppler  receivers  of  the  Institute  of  Radio  Astronomy  (Kharkov, 
Ukraine)  and  Lowell  Digisondes. 

The  aim  of  the  current  Project  investigation  is  to  develop  a  general  theory  for  the  aspect- 
sensitive  scattering  of  electromagnetic  waves  taking  into  account  ionospheric  refraction  effects  on 
the  incident  and  scattered  wave  trajectories.  We  will  specifically  analyze  a  possibility  of  signal 
channeling  in  a  given  direction  as  a  result  of  such  scattering.  The  work  on  HF  signal  trapping  in 
ionospheric  waveguides  due  to  scattering  was  pioneered  in  1970’s  (e.g.,  Erukhimov  et  al.,  1975 
Gurevich  et  al.,  1975).  Our  approach,  however,  is  quite  different  from  those  studies  as  we  will 
analyze  the  range  of  incidence  angles  at  the  lower  boundary  of  the  ionosphere  responsible  for 
channeling  signals  in  a  given  direction  rather  than  evaluate  a  “trapping  coefficient”  for  HF  waves. 
The  gyrotropic  effects  of  the  ionospheric  plasma  on  the  radio  wave  propagation  will  be  assumed 
negligibly  small.  Note  that  this  is  quite  legitimate  first-step  approach  used  in  earlier  works  as  well 
(e.g.,  Gurevich  et  al.,  1975). 

Let  a  plane  monochromatic  electromagnetic  wave  be  propagating  in  a  horizontally 
stratified  ionosphere  from  the  lower  half-space  through  the  ionospheric  layer  (see  Figure  3.1).  The 

wave  is  characterized  by  a  frequency  a  and  wave  vector  k^‘^(0)  =  {A:^‘^(0),A:|.‘^(0),k^''(0)}.  Here 
k|'^(0)  =  koSin6()Cos^Oo,  kj;'^(0)  =  kg  sin6()  sin^o  and  k^'^(O)  =  ko  cos6() ,  where  kf^=o)/c,  and 
angles  6^  and  (p^  determine  the  direction  of  the  wave  vector  at  the  lower  boundary  of  the 
ionosphere  z  =  0 .  The  wave  frequency  co  is  assumed  to  be  much  higher  than  the  electron 
gyrofrequency  cOfj ,  i.e.  co»  This  makes  it  possible  to  treat  the  ionosphere  as  an  isotropic 
medium,  and  to  assume  the  propagating  waves  to  be  transversal.  In  order  to  avoid  any  confusion, 
the  following  nomenclature  will  be  used  in  our  analysis.  The  HF  waves  will  be  regarded  as 
“incidenf’  and  “scattered”  with  respect  to  the  scattering  on  the  ionospheric  irregularities.  With 
respect  to  the  propagation  in  the  ionosphere,  the  waves  will  be  considered  as  “direcf’  and 
“reflected”. 


Distribution  A:  Approved  for  pubiic  releas;  distribution  is  unlimited. 


4/5/2015 


STCU 


PARTNER  PROJECT  P-524  EINAL  REPORT 


PAGES 


Fig.  3. 1.  Coordinate  system  for  the  problem.  X  -axis  points  toward  the  geomagnetie  north  pole. 
\  veetor  indieates  the  direction  of  the  local  magnetic  field,  (0)  is  the  incident  wave  vector  at 
the  lower  ionospheric  boundary. 


The  ionosphere  will  be  specified  as  a  plane-stratified  collisionless  dielectric  medium  containing 

random  irregularities: 

j£o(z)  +  S£(r,z),  z>0; 

£ir,z)  =  \ 

[  1,  z  <  0. 

Here  r  =  {x, y}  are  horizontal  coordinates;  £(,(z)  is  a  regular  (i.e.,  without  irregularities) 
component  of  the  dielectric  permittivity  of  the  ionosphere, 


4m^N,(z) 


where  (z)  is  a  regular  electron  density  profile,  and  e  and  m  are  the  electron  charge  and  mass. 


respectively;  and  Ss{r,  z)  = 


z) 


is  a  random  addition  due  to  electron  density  fluctuations 


SN (r ,  z) ,  characterized  by  zero  mean  <  SN (r ,  z)  >=  0  and  variance  <  [<SV (r ,  z)]^  >=  <yl  (z)  (the 
angular  brackets  <...>  stand  for  statistical  averaging).  Note  that  the  fluctuations  ds{r,z)  are 
related  to  the  relative  electron  density  variations  (W (r ,  z)  /  Nq  (z)  as: 

S£(r,  z)  =  -[1  -  (z)]SN  (F,  z)  /  No  (2)  •  (2) 

Electron  density  irregularities  SN{r,z)  (or  irregularities  in  the  dielectric  permittivity<5£’(F,z))  of 
the  magnetized  plasma  in  the  upper  ionosphere  are  highly  anisotropic,  stretched  along  the 
geomagnetic  field  direction  .  The  position  of  the  unit  vector  in  the  plane  of  the  geomagnetic 
meridian  y  =  0  (see  Figure  1)  is  specified  by  the  inclination  angle  1  (-90°  <  /  <  90°),  counted 
from  the  horizontal  plane  [Akasofu  and  Chapman,  1972].  Positive  values  of  1  correspond  to 
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pointing  downward  (northern  hemisphere),  while  the  negative  ones  correspond  to  the  upward 
direction  (southern  hemisphere).  Then,  for  the  selected  system  of  coordinates  (the  x-axis  in 
Figure  3.  1  points  to  geomagnetic  north)  we  have 

^0  =  ^oJ=  {cos/,0, -sin/}.  (3) 

In  order  to  calculate  characteristics  of  scattering  of  the  incident  electromagnetic  wave  by 
such  ionospheric  irregularities,  let  us  apply  the  single  scattering  approximation  method  [e.g., 

Rytov  et  ah,  1987].  In  this  case  it  is  assumed  that  •^<  [&•(?,  zyf  >  «  ,  and  therefore,  it  is 

possible  to  ignore  the  effects  of  the  ionospheric  irregularities  on  the  wave  trajectories  and 
calculate  them  within  the  geometrical  optics  approximation. 

Trajectory  parameters  of  the  aspect-sensitive  scattered  signals.  It  is  known  [Gershman  et  ah, 
1984]  that  HF  waves  incident  on  the  field-aligned  irregularities  of  the  upper  ionosphere  are 
predominantly  scattered  in  the  direction  determined  by  the  so-called  aspect  condition: 

(4) 

Here  K{z,)  is  the  scattering  vector,  which  is  equal  to  the  difference  of  the  wave  vectors  of  the 
scattered  k^'\z,)  and  incident  ^^‘^(z^)  fields  at  the  scattering  point  : 

=  (5) 

As  is  evident  from  Figure  3.2,  equation  (4)  is  satisfied  when  the  angles  made  by  the  magnetic  field 

vector  with  the  incident  and  scattered  wave  vectors  are  equal.  As  a  result,  the  set  of  vectors 
forms  a  characteristic  angular  cone  around  the  magnetic  field  with  the  apex  angle  v  . 


Fig.  3.2.  Scattering  geometry.  Vectors  and  are,  respectively,  the  incident  and  scattered 
wave  vectors  at  the  scattering  point  z, ;  hg  indicates  the  direction  of  the  local  magnetic  field. 

In  the  coordinate  system  shown  in  Figure  3.1,  and  with  account  of  (5),  equation  (4)  can  be 
written  as: 

■  (sin^*^^^  cos^*^^^  -  sin^^'^  cos^i9*^'’ )  +  '  (cos^^^^  -  cosd^‘^ )  =  0 »  (6) 

where  0^'^ ,  and  0^'^ ,  are  the  angles  specifying  the  wave  vector  positions  of  the  scattered 
(superscript  5 )  and  incident  (superscript  i )  fields  at  the  scattering  point  z^  (hereafter,  symbol  z^ 

is  omitted  to  simplify  the  equations).  Solving  equation  (6)  for  sin^^^^  and  cos^^*\  it  is  possible  to 
derive  the  direction  of  the  wave  vector  of  the  scattered  field.  Using  the  relation 
sin^  0^'^  +  cos^  0^^''  =  1 ,  two  possible  solutions  are  obtained: 
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COs6'i<2^  =  ,  2 


1 


^Oz+^O.COS  (p^‘ 


Xh(f^  sin6*^'^  cos^g  cos6*^'^)±/jo^  : 


(7a) 


X  cos^^'^  cos^Q  +hl^{co^  -cos^^o) 

1 


sin6»i<2^  =  71 — -2 - rm  ^i^xiKx  cos^o  +\^  cos6»^'0cos^^"^  + 

^oA^o.cos  (p^ 


7b) 


+  ^(/ig^  cos6'^‘^  cos^c>o  sin^^'^)^  +hl^{co^  cp'-'^  -cos^  cp^) 

Here  the  top  sign  stands  for  $1'^ ,  while  the  lower  sign  stands  for  .  It  is  also  taken  into  account 
that  =  (Pq  according  to  the  equations  of  geometrical  optics  for  a  plane-stratified  medium  [e.g., 
Kravtsov  and  Orlov,  1990].  Note  that  the  incident  wave  can  be  scattered  either  upward 
(cos^/2  >b),  or  downward  (cos^j^2^  <  0 ),  i.e.,  e  [0,  ;r] .  Therefore,  the  solution  space  of 

equation  (7)  with  respect  to  (p^'^  is  limited  by  the  condition: 

sin^i‘^>0.  (8) 

Another  necessary  condition  for  the  existence  of  a  real  solution  of  equation  (7)  is  a  non-negative 
argument  of  the  square  root  in  equation  (7),  i.e., 


COS^  (p^'^  >  COS^ 


cos6>^'^  cos^g  -^^sin^ 


h, 


-|2 


3(0 


h 


Ox 


(9) 


The  value  of  sin^*'^  at  is  related  to  sin^g  at  z  =  0  through  Snell’s  law  [Kravtsov  and  Orlov, 
1990] 

sin  6»g  =  Jso(z)  sin  (z)  =  Jsgizj  sin 6'^'\zJ .  (10) 

Whence  it  follows  that  if  sin^g  =  yjs^iz)  at  a  certain  altitude  z  ,  then  the  value  of  sin^*'^  is  equal 


to  1  at  that  altitude  (i.e.,  0^‘\z)  =  71/ 2).  Using  equation  (1)  this  condition  can  also  be  written  as: 

•  cos^o  =^u^(z).  (11) 

Thus,  equation  (11)  determines  the  critical  reflection  height  z^^  for  the  ray  with  the  incidence 
angle  6^  at  the  lower  boundary  of  the  ionospheric  layer.  This  is  why  when  we  determine 
cos0^'\z^)  two  situations  are  possible,  depending  on  the  ratio  between  rycos^g  and  (see 
Figure  3.3).  If  <ycos^g  >  <»^^ ,  then  the  wave  penetrates  through  the  ionospheric  layer  (no 

reflection).  In  this  case,  within  the  entire  volume  of  the  ionospheric  layer  with  irregularities,  the 
scattering  occurs  only  at  the  ascending  part  of  the  wave  trajectory  (the  direct  wave  is  scattered 
alone),  and  the  value  of  cos^^'^(zJ  according  to  [Kravtsov  and  Orlov  (1990)]  is: 

cos(9^'^(z,)  =  ^/cos^"^J^^a3^(zj7^.  (12) 

If  rucos^Q  <  ,  then  the  wave  is  reflected  at  z  =  z^, ,  determined  by  equation  (1 1),  and  the  range 

of  z,  is  limited  as  0  <  z,  <  z^^ .  As  a  result,  two  waves  arrive  at  the  scattering  point,  one  being 
scattered  on  the  ascending  part  of  the  trajectory  (direct  wave),  the  other  on  the  descending  path 
(reflected  wave).  Then  for  cos0^'\zX  g^ts: 

cos6>^''^ (z,)  =  ±^ cos^  6>o  - (ol (z, )/co\  (13) 

where  “+”  and  stand  for  the  direct  and  reflected  waves,  respectively. 

1 . 
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a)  b) 

Fig.  3.3.  Trajectories  of  the  incident  wave  for  (ocos6q>  {6^=6^^)  and  cocosO^kco^^ 
(^0  =  ^02) »  where  co^  is  the  critical  frequency  of  the  layer  (plasma  frequency  at  the  height  of  the 
layer  maximum  z„).  In  the  first  case  (a),  the  direct  wave  =  alone  reaches  the 

scattering  height  ,  while  in  the  second  case  (b)  two  waves  reach  the  height  z^ ,  specifically,  the 
direct  wave  ( 6^'^  (z^ )  =  9^1 )  and  the  wave  reflected  at  z„  ( 9^'^  (z^ )  =  9^^2  )• 

Now,  consider  the  scattered  wave.  Let  the  background  electron  density  profile  Nq{z)  be  a  smooth 
function  with  a  single  maximum  N^{z^)  at  z„  (single  layer  model).  If  <x)^^><x) ,  then  the  waves 
scattered  downward  ( cos,9^'^  <  0 )  and  upward  ( cos^^*^  >  0 )  both  reach  the  lower  boundary  of  the 
ionospheric  layer  at  z  =  0  with  an  angle  0^^  "'  G[2r/2,;r],  which  can  be  found  from  an  equation 
similar  to  equation  (10),  viz. 

AneZ’  =  .  (14) 

Ifm,,  <  CO ,  then  only  the  waves  scattered  at  angles  ^^^^(z^)  >  9*^^  can  reach  the  lower  boundary  of 
the  ionosphere  (z  =  0).  Here  9^^^  e[0,;T/2]  is  a  certain  critical  angle  whose  value,  with  account 
of  equation  (1)  and  Snell’s  law  (10),  can  be  determined  as: 


9li^ ,  the  scattered  wave  goes  into  the  upper  half-space  and  then  penetrates  through  the 

ionosphere.  Note  that  in  the  case  of  a  multi-layer  ionosphere  (e.g.,  E  -  and  F  -layers  are  present) 
the  waves  scattered  in  certain  directions  can  be  “trapped”  by  the  interlayer  ionospheric  waveguide. 
Such  a  possibility  will  be  investigated  in  the  next  section,  but  first  let  us  analyze  the  amplitude 
characteristics  of  the  aspect-sensitive  scattered  signals.  Scattering  cross-section.  The  effective 
differential  cross-section  of  a  random  medium  is  used  to  characterize  the  energy  scattered  by  a 
unit  volume  into  a  unit  solid  angle  in  a  given  direction  with  a  unit  flux  density  of  the  incident 
radiation  [e.g.,  Rytov  et  ah,  1987].  In  the  approximation  of  a  single  scattering  of  electromagnetic 
waves  in  random  isotropic  medium  with  a  regular  dielectric  permittivity  =  const,  the  scattering 
cross  section  is  [Rytov  et  ah,  1987;  Gershman  et  ah,  1984]: 
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(16) 

Here  is  a  three-dimensional  spatial  spectrum  of  fluctuations  de',  K{z^)  is  the 


Q{K{z,))  =  P{z,)^<^,{K{z,))  . 


scattering  vector  given  by  equation  (5);  and  P{z,)  =  1- 


is  a  polarization  factor, 


where  e^'\z^)  is  the  unit  vector  of  the  electric  field  polarization  of  the  incident  wave.  Using  the 

same  geometrical  optics  assumptions  as  in  Rytov  et  al.  [1987],  it  can  be  shown  that  equation  (16) 
is  applicable  to  the  analysis  of  the  scattering  of  a  single  quasi-plane  wave  in  the  medium  with 
smooth  spatial  variations  in  £'o(Zj).  The  factor  Piz,)  is  dependent  on  the  incident  wave 
polarization.  For  instance,  according  to  [Rytov  et  ah,  1987]  in  the  case  of  a  linear  polarization 

P(z,)  =  sin^;if(zJ,  (17) 

and  for  a  circular  polarization 


R(z,)  =  ^(l  +  cos^  riz,)). 


(18) 


where  xiz^)  /(z,)  denote  the  angles  made  by  the  scattered  field  wave  vector  k^^\z,)  with 

the  polarization  vector  e^’\z^)  and  wave  vector  k^‘\z,)  of  the  incident  field,  respectively. 

The  anisotropic  power  law  model  (e.g.,  Gershman  et  al.  [1984])  is  a  conventionally  used 
approximation  for  the  spatial  spectrum  O^(^)  of  the  fluctuations  S£(r,z)  in  the  magnetized 
plasma  of  the  upper  ionosphere  for  the  inertial  interval  of  the  wave  numbers.  In  this  model, 
0^(^(zj)  can  be  written  as 


(19) 


where  C^(z)~cr^(z)  is  a  normalization  factor,  with  cr^(z)  representing  the  variance  of  the 
fluctuations  Ss;  and  are  longitudinal  and  transversal  (with  respect  to  h^)  components  of 

the  scattering  vector  ^  and  are  characteristic  longitudinal  and  transversal  external  scale- 
sizes  of  the  ionospheric  turbulence  with  Lj_;  and  3  <  p  <  4 . 

If  condition  (4)  is  satisfied,  then  equation  (16)  yields 

1  ^  ^  L  ^  ^  l-p/2 


Q(Zs)  =  P(Zs )  ^  CAz,  )[l  +  Kl  (z,  )L\  ] 


where 


Kl{z,)  =  kl 


sin^  Oq'*  -i-sin^  0^^  -i-(sin0o''’  cos^p**^  -sin^g  cosq)^'''f 


Oz 


—  2sin0o''*  sin^o  cos(^*'* 


(20) 


(21) 


Equation  (21)  is  inapplicable  in  the  specific  case  of  the  horizontal  magnetic  field,  =0,  (at  the 
magnetic  equator).  It  can  be  shown  that  in  this  case  K]_{z^)  becomes 

=  -l-sin^  -2  sin  sin^o  cos(^*'^ 

Using  equation  (2),  C^(z^)  can  be  expressed  as 


Distribution  A:  Approved  for  pubiic  releas;  distribution  is  unlimited. 


4/5/2015 


STCU 


PARTNER  PROJECT  P-524  PINAL  REPORT 


PAGE  13 


Cl  (z,  )~C7l(z,)  =  (l-£o  (Zs  )y 


®p(z,)  cjI(z,) 


(22) 


Nliz,)  iVo'(^.)  ■ 

Therefore,  from  equation  (20)  it  is  possible  to  determine  the  scattering  cross-section  in  the 
direction  defined  by  the  aspect  condition  (4)  for  a  plane  electromagnetic  wave  propagating  in  an 
ionospheric  layer  with  anisotropic  electron  density  irregularities.  Note  that  since  (p^'^{z)  =  (p^  and 

(z)  =  (Pq^  ,  all  the  terms  in  equation  (20)  can  be  expressed  through  the  position  angles  of  the 


incident  and  scattered  field  wave  vectors  at  the  lower  boundary  z  =  0  of  the  ionospheric  layer. 
Discussion.  As  an  example,  let  us  first  consider  scattering  of  a  plane  circularly  polarized 
electromagnetic  wave  propagating  vertically  in  an  ionospheric  layer  with  anisotropic  irregularities. 
The  wave  frequency  is  assumed  to  be  lower  than  the  critical  frequency,  co<co^^.  Note  that  these 


are  typical  conditions  for  the  HF  heating  of  the  ionosphere. 

In  a  single  layer  model  ionosphere  two  components  will  be  scattered  in  the  height  range 
0<Zj  <z„  {z„  is  determined  from  (11)  with  cos6>q  =1),  which  are  the  direct  and  reflected 
waves.  Then  equation  (7)  with  account  of  equation  (13)  yields: 


cosOlfl 


^hi-hlcol^ 
hi,  +  hi  cos' ’ 


(23a) 


sinC=  + 


2hoAz  ^os^^'^ 


hlz 


+  h 


Ox 


cos'  (p^'^ 


>0, 


(23b) 


where  the  “+”  and  signs  stand  for  the  scattering  components  of  the  direct  (subscript  d )  and 
reflected  (subscript  r )  waves,  respectively. 

As  can  be  seen  from  equation  (23b),  in  the  northern  hemisphere  (  <  ^ )  the  direct  wave  is 


scattered  southward  ( -  90°  <  <  90° ),  while  the  reflected  is  scattered  northward 


(90°  <  <  270°) .  In  the  southern  hemisphere  the  situation  is  reversed.  Near  the  magnetic 

equator  (/j^  «0),  as  follows  from  equation  (7b),  both  components  are  scattered  predominantly 
within  the  plane  perpendicular  to  the  geomagnetic  field  lines  at  all  angles  e  [0,  k}  . 


To  analyze  the  orientation  of  the  scattered  field  wave  vector  at  the  exit  from  the 
ionospheric  layer  it  is  convenient  to  introduce  an  angle  =71-0^  ^ .  Then,  according  to  Snell’s 
law  and  using  equations  (23b)  and  (1),  one  gets 


sin^(^> 


2,/r 


(24) 


As  an  example,  Figure  4  shows  isolines  for  (p^'^  =  75°(255°)  ,45°(225°)  and  0°(180°), 


which  are  plotted  in  the  coordinates  a  =  G)l{z^')l (o'  and  ■  The  inclination  angle  I  of  the 


geomagnetic  field  is  set  equal  to  77.5°  (recall  that  /io^=cos/  and  /iQ^=-sin7),  which 

corresponds  to  the  location  of  the  EISCAT  HF  heater  (69°35'N,  19°14'E).  In  this  case,  the  wave 
vector  of  the  scattered  component  at  the  lower  boundary  of  the  ionosphere  ( z  =  0 )  can  be 
significantly  off  the  vertical  despite  strongly  vertical  propagation  of  the  primary  wave.  The  value 
of  this  offset  increases  as  a  decreases,  (i.e.,  with  decrease  of  the  scattering  altitude  z, ),  and/or  as 
the  scattering  direction  approaches  the  plane  of  the  magnetic  meridian.  In  this  example  the 
maximum  offset  is  equal  to  71  —  21  =  25° ,  for  z^  =  0  and  (p^'^  =  0,  .  Analysis  of  equation  (24) 

shows  that  such  a  dependence  is  valid  for  the  geomagnetic  field  inclinations  |  /  |>  /  4 .  Otherwise, 

if  I  /  |<  .TT  /  4 ,  the  maximum  offset  of  the  aspect-sensitive  scattered  wave  from  the  vertical, 
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(i.e.,  the  scattering  within  the  horizontal  plane)  will  occur  for  and 

=±^0z/V- 


Fig.  3.  5.  Angular  distributions  of  the  value  \Q\og\Q{z,) ! calculated  using  equation  (25) 
with  p  =  1 1/3  for  L^l \=\  (a)  and  LJX,=5  (b).  The  solid,  dashed  and  dotted  lines 
correspond  to  the  levels  a  =  co^p{z,) ! co^  =0.3,  a  =  0.1 5  and  «  =  0.99,  respectively.  Calculations 

were  made  for  the  vertical  incidence  at  the  latitudes  corresponding  to  the  EISCAT  location 
(7  =  77.5°).  The  value  of  the  scattering  cross-section  increases  with  a  (i.e.,  with  the  scattering 

height  ),  reaching  its  maximum  at  a  =  1 ,  which  corresponds  to  =  0° . 


To  calculate  the  scattering  cross  section  in  the  direction  defined  by  equations  (23)  we  combine 
equations  (18),  (21),  and  (23)  and  Snell’s  law  (equation  (10))  obtaining  for  the  vertically 
propagating  wave 


P(z,)- 


K,  +Cg°sV*' 

^  +hl..cos^ 


(h, 


ilcos'  r 

:.2 _ 2 


Kl(zY  = 


4^ogo(z,)/?o,cos  (p' 

hL  +  hL  cos^ 


/•') 


‘Oz  '  “-Ox 

Substitution  of  these  expressions  together  with  equation  (22)  into  (20)  yields 

.4/_  ^  ,_4  ,  ,.4 _ 4  (s) 


Q(zY~Qo(Zs)- 


<^p(z,)  +hQ^cos^  <p 


CO 


ih 


^0z+ Ceos'  cp^^^f 


1+ 


d^o^oAgpC^Jcos  cp' 

C+Ccos'  A 


G) 


-pH 


(25) 


where  Qq{z,)  = 


tAq  <7j^ 


o 


2  NlizY 


.  For  quiet  ionospheric  conditions  it  is  usually  assumed  that  the 


plasma  density  fluctuations  are  proportional  to  the  background  electron  density,  i.e.. 


Nl(z,) 


=  const  =  A  •  Thus  one  can  write  2o  (z^ )  = 


nk, 


0_^2 
^jv  ■ 
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Figure  3.5  shows  the  angular  distribution  of  the  value  101og[2(Zj)/2o(Zj)],  calculated  with 
equation  (25)  with  p  =  \\l3  and  Lj^l A^=\  (top  panel),  and  LJ\=5  (bottom  panel)  for  the 
vertical  incidence  case  for  the  latitude  corresponding  to  the  EISCAT  location  (  /  =  77.5°).  Here 
Xq  =  f !  c  =  271 !  is  the  free  space  wavelength.  The  results  are  shown  for  the  scattering  of  the 

reflected  wave  (scattering  in  the  northward  direction,  e  ±90° ).  For  the  direct  wave  (scattering 
in  the  southward  direction,  el 80° ±90°)  the  angular  distributions  are  completely  symmetric. 
In  order  to  be  able  to  at  least  roughly  relate  these  data  to  the  scattering  heights  (or  the 
corresponding  plasma  frequency  co^iz,))  Figure  5  shows  isolines  for  a  =  co^p{z,) ! co^  =  0.3  (solid 
line),  a  =  0.7 5  (dashed  line),  and  cir  =  0.99  (dotted  line).  Note  that  and  oj^iz,)  are 

related  through  equation  (24).  As  is  evident  from  Figure  5,  the  value  of  the  scattering  cross-section 
increases  with  a  (i.e.,  with  the  scattering  height  z,),  reaching  its  maximum  at  a  =  l,  which 

corresponds  to  =0.  This  is  because  cr^(z^)  increases  with  height  (since  it  is  assumed  that 
crj/iz,)  ~  Nq(z^)  ~  cOp(z^)),  and  also  because  the  second  term  in  the  square  brackets  of  equation 
(25)  gets  smaller  (since  50(25)  ^  0  as  ^  z„ ).  The  latter  effect  can  be  treated  as  an  increase  of 
the  signal  wavelength  in  the  ionosphere  as  it  approaches  the  reflection  point  since 
k{z,)  =  271!  7l{zJ  =  •  Apparently,  this  is  also  the  reason  for  the  widening  of  the 

azimuthal  distribution  of  the  scattering  cross-section  with  the  decrease  of  .  For  example, 
already  for  *5°...10°  the  scattering  cross-section  is  practically  isotropic  in  the  azimuthal 
plane,  and  smaller  values  lead  to  greater  values  of  the  scattering  cross-section.  This  would 

suggest  that  in  the  presence  of  several  ionospheric  layers  (e.g.,  E  and  F)  a  certain  fraction  of  the 
scattered  wave  energy  can  be  captured  by  the  interlayer  ionospheric  waveguide.  Since  such 
waveguide  propagation  is  characterized  by  a  very  small  attenuation,  it  is  quite  plausible  that  this 
mechanism  supported  the  super  long  distance  propagation  of  the  EISCAT  signals  observed  in  the 
“self-scattering”  experiment  [Zalizovski  et  al.,  2009].  The  potential  of  exciting  the  ionospheric 
interlayer  waveguide  by  the  aspect-sensitive  scattered  HF  signals  requires  greater  investigation.  It 
should  be  noted  that  the  role  of  scattering  in  HF  radio  wave  trapping  into  ionosphere  waveguides 
was  rather  intensively  examined  as  early  as  in  1970’s  [e.g.,  Gurevich  et  al.,  1975;  Erukhimov  et 
al.,  1975].  However,  in  contrast  to  the  previous  studies,  we  will  be  primarily  interested  in 
determining  the  range  of  0^  and  (p^  responsible  for  channeling  scattered  signals  in  a  given 

direction  rather  than  in  estimating  the  trapping  coefficient  for  an  incident  wave. 

Let  the  ionosphere  be  represented  by  two  layers,  e.g.,  E  and  F ,  each  characterized  by 
their  respective  critical  frequencies  and  co^^p  and  critical  heights  z^p  and  z„f  .  We  denote 

the  minimum  plasma  frequency  inside  the  E-F  valley  at  the  height  z„  as  co^ .  Further  assume  that 
co^^P  >  CO,  i.e.,  the  wave  with  frequency  co  is  reflected  from  the  F  -layer  (the  upper  “wall”  of  the 
waveguide)  at  any  incidence  angle.  Obviously,  the  interlayer  waveguide  can  be  excited  by  the 
waves  scattered  in  the  height  range  z„p  <  z,  <  z„p  only.  Therefore,  it  is  necessary  for  the  direct 
wave  to  penetrate  the  E  -layer,  i.e., 

cocos0g  >  co^p .  (26) 

As  follows  from  Snell’s  law,  in  order  for  the  scattering  component  to  be  confined  inside 
the  waveguide,  i.e.,  to  be  reflected  from  the  E-layer  (the  lower  “wall”  of  the  waveguide),  the 
following  condition  must  be  met: 

^^o(Zs)sm0^'\zJ>  V^o(Zcr£)  (27a) 
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which,  with  account  of  equation  (1),  can  be  recast  in  a  form  similar  to  equation  (15): 

sin6^^"^(zj  > 


CO  -CO, 


crE 


(27b) 


where  is  determined  by  equation  (7b).  Since  sin^^^^(z^)  <  1 ,  it  follows  from  equation 

(27b)  that  the  scattering  should  occur  within  the  height  range  <  z,  <Zp  ,  where  Zp  e  [z^ ,  z^p  ] , 


and  can  be  determined  from  cOp(zp)  =  co^^p  ■ 

Further  analysis  will  be  carried  out  for  the  specific  conditions  of  the  “self-scattering” 
experiment  conducted  by  Zalizovski  et  al.  [2009].  We  are  primarily  interested  in  the  possibility  of 
channeling  the  EISCAT  signals  inside  the  interlayer  waveguide  in  the  directions  of  the  receive 
sites  located  at  UAS,  SPB,  and  RAO.  In  the  chosen  coordinate  system,  «135°  for  the  UAS 
station  and  ~  225°  for  the  other  two  stations.  In  most  heating  sessions  the  EISCAT  heater 
transmitted  toward  magnetic  zenith  at  one  or  two  close  frequencies  around  4.04  MHz.  The 
measurements  from  the  collocated  ionosonde  showed  the  presence  of  strong  E  layers,  with  the 
critical  frequencies  close  to  and  even  above  the  heating  frequencies.  For  the  analysis 
fcrE  =3.9  MHz  and  « 1.4 MHz  are  assumed  and  the  calculations  are  not  limited  to  the  vertical 

incidence  case,  making  it  possible  to  investigate  the  effect  of  the  angles  of  direct  wave  incidence. 
As  known  [Gurevich,  2007],  artificial  ionospheric  turbulence  is  effectively  produced  by  the  HF 
heating  in  the  height  range  with  the  lower  boundary  Zp^ ,  determined  from  the  upper  hybrid 

resonance  condition  =  fp(ZpH)  +  fniZun)  ’  where  is  the  gyrofrequency;  and  the  upper 
boundary  z„  ,  determined  by  the  wave  critical  reflection  condition  equation  (11).  For  this  reason 
the  plasma  frequency  fp(z^)  in  the  calculations  is  assumed  to  vary  within  the  scattering  region 

from  fp(ZpH)  =  ^f^-fH(ZpH)  « 3.79MHz,  to  fp(z,,)  =  f,,E  =3.9 MHz.  Figure  6  shows  results 
of  the  calculations  for  such  conditions  in  the  system  of  coordinates  with  6^  (counted  radially),  and 
(counted  counterclockwise).  The  plots  show  the  distributions  of  sin^^^*  values  which  satisfy 
the  conditions  of  excitation  of  the  interlayer  ionospheric  waveguide  (see  equation  (27b)).  The 
calculations  were  made  for  /■  =  4.04MHz,  ^^*^=135°  and  /),=3.79MHz  (Figure  6a)  and 
fp  =  3.87MHz  (Figure  6c).  Figures  6(b)  and  6(d)  show  the  respective  distributions  of  the  value 
101og[2(zJ/2o(zJ]  for  p  =  ll/3  and  L_^l \=\.  The  values  are  given  in  grayscale  (see  the 

legends  on  the  right).  As  follows  from  equations  (7)  and  (21),  the  distributions  of  sin^^^'  for 
=  225°  will  be  the  same  as  in  Figure  6  (a)  and  (c),  while  the  distributions  of 
101og[2(zJ/2o(zJ]  will  look  like  mirror  images  of  Figure  6  (b)  and  Figure  6  (d)  with  the 
reflection  line  at  =180°.  Note  that  limits  of  the  region  from  which  the  scattered  signals  are 

JC^2  ^2 

— ^ (see 

equation  (27b)),  which  is  approximately  equal  to  0.754  for  =  3.19  MHz  and  0.909  for 
fp  =3.87  MHz.  The  dashed  line  in  Figure  6  shows  the  EISCAT  antenna  beam  at  the  half-power 
level.  As  can  be  seen,  the  signal  propagation  through  the  interlayer  duct  channel  in  the  azimuthal 
directions  =135°  and  cp^^^  =225°  is  possible  for  the  given  orientation  of  the  heater  antenna 
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only  for  a  very  narrow  range  of  angles  near  6^  « 1 5°  and  ^(,=180°,  with  characteristic  sizes  (for 
fp  =  3.87  MHz)  of  ±15°  in  azimuth  and  less  than  \°  m  0^  (see  Figure  3.6  a). 


a) 


b) 


c)  d) 

Fig.  3.  6.  Distributions  of  the  sin6>^*^  values  which  satisfy  the  condition  of  excitation  of  the 
ionospheric  interlayer  waveguide  (equation  (27b)).  Calculations  were  made  for  f  =  4.04  MHz, 

=135°  and  =  3.79  MHz  (a)  and  =  3.87 MHz  (c).  Panels  (b)  and  (d)  show  the  respective 
distributions  of  the  values  101og[2(zJ/2Q(zJ]  for  p  =  ll/3  and  L^IXq=\.  The  values  are 
shown  in  grayscale  (see  legend  on  the  right  side)  in  the  coordinate  system  with  6^  (shown  as 
radius)  and  %  (counted  counterclockwise).  The  dashed  line  shows  the  EISCAT  antenna  beam  at 

the  half  power  level.  It  can  be  seen  that  for  the  given  orientation  of  the  heater  antenna,  signal 
propagation  in  the  interlayer  waveguide  toward  the  receive  sites  is  possible  only  for  a  very  narrow 
range  of  angles  near  » 15°  and  =180° ,  with  characteristic  sizes  (for  =  3.87  MHz)  about 

±15°  in  azimuth  and  less  than  1°  in  . 

Accordingly,  the  linear  size  of  the  scattering  region  at  200  km  is  approximately  45  to 
50  km  across  and  3  to  4  km  along  the  magnetic  meridian.  Note  that,  as  follows  from  equation  (7), 
this  case  corresponds  to  the  scattering  of  the  ionospherically  reflected  signal.  The  sharp  cutoff  at 
6q  » 1 5°  is  determined  by  condition  (26)  for  the  direct  wave  to  penetrate  the  E  -layer.  The  value 

of  101og[g(zJ/2o(2j]  in  this  angular  range  varies  insignificantly,  from  about  -4  to  -12  (see 
Figure  3.6  b).  As  increases,  i.e.,  as  the  scattering  height  approaches  the  critical  reflection 
height  z„  ,  the  angular  and  linear  sizes  of  the  scattering  region  decrease  by  approximately  a  factor 

of  2  (see  Figure  3  6  c),  while  the  cross-section  changes  insignificantly  (see  Figure  3.  6  d). 

The  above  analysis  supports  the  mechanism  of  the  “self-scattering  effect”  that  was 
suggested  by  Zalizovsky  et  al.  [2009]  for  the  explanation  of  their  experimental  results  of 
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monitoring  the  EISCAT  transmission  at  three  remote  locations  (SPB,  RAO,  and  UAS).  Recall  that 
during  the  experiment  the  variations  of  the  Doppler  frequency  shifts  recorded  at  all  the  sites  were 
almost  identical.  During  the  experiment  significant  variations  of  the  electron  density  were  also 
observed  by  the  EISCAT  incoherent  scatter  radar  in  the  height  range  from  100  to  200  km.  Such 
variations  can  significantly  affect  the  Doppler  frequency  shift  and  amplitude  of  the  HE  heating 
signal  as  was  demonstrated  by  the  authors  using  computer  modeling.  With  the  above  analysis  we 
have  demonstrated  that  the  angular  sizes  of  the  scattering  region  responsible  for  the  excitation  and 
subsequent  channeling  of  the  HE  signals  through  the  interlayer  waveguide  to  the  receive  sites  are 
rather  small  (see  Eigure  3.6).  Therefore,  it  is  reasonable  to  assume  that  all  remotely  recorded 
signals  first  travel  through  the  same  irregularities  in  the  lower  ionosphere,  and  are  equally  affected 
by  them.  Then,  the  signals  reflected  from  the  ionospheric  E-layer  are  scattered  by  field-aligned 
irregularities  produced  by  the  same  EISCAT  heating  transmission  and  propagate  toward  the 
receive  sites  through  the  ionospheric  waveguide.  Since  the  scattering  region  guiding  the  signals 
into  the  interlayer  duct  charmel  for  propagation  into  different  directions  is  the  same,  it  is  plausible 
to  expect  that  the  variations  of  the  Doppler  frequency  shifts  and  signal  amplitudes  recorded  at 
greatly  dispersed  receive  sites  will  be  very  similar,  as  was  observed  during  the  “self-scattering” 
experiment  reported  by  Zalizovski  et  al.  [2009]. 

Conclusions,  hr  this  chapter  the  ray  optics  approximation  and  small  perturbation  method  were 
used  to  analyze  the  effect  of  refraction  on  the  scattering  characteristics  of  HE  signals  scattered  by 
random  field-aligned  irregularities  of  the  ionospheric  plasma.  The  calculations  are  carried  out  for 
an  isotropic  plane-stratified  (on  average)  ionosphere,  i.e.,  the  incident  and  scattered  waves  were 
both  assumed  to  propagate  along  unperturbed  trajectories.  The  equation  of  the  so-called  aspect- 
sensitive  scattering  is  derived  which  relates  the  trajectory  characteristics  of  the  incident  and 
aspect-sensitive  scattered  waves.  The  scattering  cross-section  is  calculated  within  the  Bom 
approximation  for  a  plane  electromagnetic  wave  propagating  in  the  ionospheric  layer  with 
anisotropic  electron  density  irregularities.  Eor  the  vertical  incidence  scenario  the  intensities  of  the 
aspect-sensitive  scattered  wave  (cross-section)  and  its  exit  angle  are  calculated  as  functions  of  the 
scattering  direction  and  height.  A  possibility  of  excitation  of  the  interlayer  ionospheric  waveguide 
by  the  aspect-sensitive  scattered  signals  with  a  selected  orientation  of  the  horizontal  projection  of 
the  wave  vector  is  investigated  in  dependence  on  the  angles  of  direct  wave  incidence.  It  is 
demonstrated  that  such  possibility  can  be  implemented  using  powerful  HE  heaters,  like  EISCAT 
or  HAARP.  Under  certain  conditions  (e.g.,  specific  orientation  of  anterma  beam),  the  scattered  HE 
signals  can  be  trapped  inside  the  waveguide  and  subsequently  channeled  to  very  long  distances 
because  of  the  very  small  attenuation  characteristic  of  the  interlayer  waveguide. 

Thus,  the  aspect-sensitive  scatter  mechanism  can  explain  the  “self-scattering  effect”  which 
was  observed  using  the  EISCAT  transmission  by  Zalizovsky  et  al.  [2009].  Therefore,  there  exists 
the  potential  for  the  practical  applications  in  HE  communication,  although  further  studies  are 
required  for  a  more  reliable  assessment.  In  particular,  it  is  necessary  to  consider  the  suggested 
mechanism  with  account  of  the  magnetic  field  effects  which  can  influence  essentially  propagation 
of  the  incident  and  scattered  waves.  Addressing  this  issue  properly  would  require  application  of 
full  scale  3D  numerical  ray-tracing.  The  authors  plan  to  carry  out  the  respective  analysis  in  the 
future  and  present  the  results  in  a  separate  publication.  As  a  first  approximation,  however,  we 
think  that  neglecting  the  geomagnetic  field  effect  on  radiowave  trajectories  is  a  legitimate 
approach,  which  was  also  taken  by  other  authors  [e.g.,  Gurevich  et  al.,  1975].  The  results  of  this 
investigation  can  be  useful  for  the  analysis  and  interpretation  of  the  experiments  on  HE  signals 
scattering  from  artificial  and  natural  ionospheric  irregularities  in  the  ionosphere,  for  the 
development  of  new  methods  of  diagnostics  of  the  near-Earth  plasma  turbulence,  and  also  for  the 
investigation  of  the  mechanisms  of  super- long  range  propagation  of  HE  electromagnetic  emission. 
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4.  Studies  of  the  Ionospheric  Turbulence  Excited  by  the  Fourth  Gyroharmonic  at 
HAARP  and  ionospheric  wave  guide  excitation 

There  are  the  results  of  experiments  conducted  during  the  HAARP  BRIOCHE  June  2014 
campaign,  whose  objective  was  to  study  the  development  of  artificial  ionospheric  turbulence. 
During  the  experiments,  the  heating  frequency  was  stepped  up  near  the  4*  gyroharmonic,  and  the 
power  of  the  heating  HE  radiation  was  varied.  Our  diagnostics  included:  measurements  of  phase- 
derived  Slant  Total  Electron  Content  using  the  L1/L2  GPS  signals  received  at  HAARP; 
measurements  of  Stimulated  Electromagnetic  Emission  (SEE)  conducted  15  km  away  from  the 
HAARP  site;  ionograms  from  HAARP ’s  digisonde  and  reflectance  data  from  Kodiak  radar;  and 
detection  of  the  HAARP  HE  radiation  scattered  into  the  ionospheric  channel  and  propagated  for 
15.6  Mm  to  the  receiver  at  Antarctica  Peninsula. 

Additionally,  a  new  method  of  diagnostic  was  introduced  during  this  campaign.  This 
method  was  developed  of  the  IRA  NASU  team  in  a  frame  of  the  EOARD  -  STCU  Partner  Project 
P-524.  HAARP  radiated  HE  waves  can  be  scattered  into  the  ionospheric  waveguide  by  artificial 
irregularities,  and  can  propagate  along  the  waveguide  like  a  sound  wave  in  the  whispering  gallery 
mode.  The  waveguide  was  oriented  along  the  Earth’s  terminator  which  was  passing  over  the 
transmitting  and  receiving  sites  simultaneously.  Furthermore  the  electron  density  gradients  formed 
in  the  ionosphere  when  the  terminator’s  passage  occurs  near  the  receiving  site  could  scatter  the  HE 
signals  from  the  waveguide  to  the  groundbased  receiver.  We  were  able  to  detect  the  scattering 
HAARP  signals  on  the  ground  at  Ukranian  Antarctic  Station  (UAS)  (coordinates  65.25  S,  64.25 
W)  at  15.6  Mm  from  HAARP.  The  HE  signals  were  recorded  within  frequency  band  500  Hz 
centered  at  the  carrier  frequency  of  the  HAARP  pumped  signal.  The  data  acquisition  system 
collected  records  of  the  intensity  and  the  Doppler  spectra  of  the  signal  with  time  resolution  1  sec 
and  5  sec.  This  paper  introduces  a  successful  attempt  to  use  the  effect  of  long  distance  HE 
propagation  for  diagnostic  of  the  artificial  turbulence.  Finally,  all  the  experiments  were  diagnosed 
by  the  iononosonde,  and  by  the  Kodiak  coherent  radar  located  670  km  South  West  from  HAARP. 
Ionospheric  conditions.  The  two  experiments  discussed  were  conducted  during  20  minutes  each 
starting  at  about  3  UT,  i.e  around  7  p.m.  local  time.  During  the  first  day  (06/06/14)  the  ionosphere 
was  slightly  disturbed  (5B~50  nT);  a  noticeable  sporadic  E-layer  existed  (foEs~4.5  MHz);  foF2  was 
in  the  frequency  range  5. 6-5. 7  MHz  while  the  F2  peak  was  located  at  hmF2=270  km.  Considering 
the  oblique  propagation  of  the  heating  wave,  its  frequency  was  reflected  from  the  F2  peak  thus,  for 
the  majority  of  the  experiment,  the  heating  wave  was  strongly  absorbed  in  the  ionosphere  due  to 
anomalous  absorption.  An  exception  was  the  time  around  3:04  UT  when  the  ionogram  revealed 
that  foF2  dropped  below  5.4  MHz,  thus  the  heating  frequency  exceeds  the  critical  frequency,  and 
the  ionosphere  became  transparent  to  the  HE  wave. 

During  the  second  day  (06/07/14)  the  ionosphere  was  also  quiet  (5B~30  nT)  with  Eg  layer 
present.  However,  the  ionosphere  was  unstable,  half  of  the  ionograms  reveal  that  foF2  fell  to  4.05- 
4.5  MHz,  thus  the  HE  frequency  was  higher  than  the  critical  frequency  and  absorption  was  low 
and  the  artificial  turbulence  was  almost  absent,  while  the  rest  of  the  ionograms  show  foF2>5.7 
MHz  and  hmF2=270km,  and  thus  at  the  time  the  artificial  turbulence  was  pumped.  Note  that  due  to 
the  experiment  restrictions  the  minimum  heating  frequency  was  5.67  MHz  which  was  above  the 
4*  gyro  resonance.  In  our  March  2013  experiment  [Najmi  et  al.,  2014]  we  have  obtained  the 
resonance  frequency  5.76  MHz  from  the  SEE  spectrum,  at  which  frequency,  the  DM  disappeared. 
The  ionogram  shows  that  5.76  MHz  waves  are  reflected  at  190  km. 
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An  estimate  based  on  the  dipole  model  of  the  geomagnetic  field  shows  that  in  the  June  2014 
experiments  the  4*  gyro  resonance  occurs  at  5.6  MHz,4/,g  =  5.16MHz 


^  R^+mkrn^^ 
+  250km 


,  where 


250  km  is  the  reflection  height  of  5.6  MHz  taken  from  ionograms.  Thus  the  lowest  heating 
frequency  5.76  MHz  is  close  to  the  BUM  cutoff  which  is  about  8  kHz  above  the  4*  gyroharmonic, 
i.e.  at  5.68  MHz  [Leyser  et  ah,  1994]. 

SEE  observations.  Stimulated  Electromagnetic  Emission  (SEE)  signals  were  measured 
simultaneously  with  the  STEC  using  an  HE  detector  operated  by  the  Naval  Research  Laboratory 
15  km  away  from  the  HAARP  site.  Figure  4.1  shows  the  power  spectral  density  (PSD)  of  the  SEE 
emission  for  each  of  four  chosen  heating  frequencies,  5.67,  5.70,  5.73  and  5.76  MHz,  which 
reveal  distinct  BUM.  At  higher  heating  frequencies  BUM  disappears.  The  traces  are  averaged  over 
a  10s  portion  of  the  heating  period  with  constant  ERP.  The  heating  frequency  is  shown  by  the 
highest  peak  in  the  center  at  Af=0,  the  down  shifted  maximum  (DM)  is  on  the  left  side  of  the 
heating  frequency  at  Af=-10  kHz,  while  the  broad  upshifted  maximum  (BUM)  is  on  the  right  side 
of  the  plots  in  the  range  of  Af=+30-130  kHz.  In  each  of  the  panels  of  Fig.  4.1,  the  blue  trace 
corresponds  to  the  power  2.5  kW,  green  to  5.0  kW  and  red  to  10  kW  per  transmitter  (the  power 
level  25,  50,  and  100%  respectively).  Note  that  the  shape  of  SEE  signal  saturates  at  the  half  full 
power.  Figure  5a,b  shows  the  amplitude  of  the  PSDs  normalized  by  its  peak  values  as  a  function 
of  heating  power.  Here  Fig.  4.2  a)  corresponds  to  the  BUM  amplitude  while  Fig.  4.2  b)  shows  the 
evolution  of  the  DM  amplitude.  In  those  two  cases  saturation  at  90%  level  occurs  at  about  ‘A  of 
the  full  power. 


PSD  over  Different  ERP  06-Jun-2014;  Pump  =  5.67  MHz 


PSD  over  Different  ERP  06-Jun-2014;  Pump  =  5.70  MHz 


Fig.  4.1.  Power  spectral  density  (PSD)  of  the  SEE  emission  for  each  of  four  chosen  heating 
frequencies,  5.67,  5.70,  5.73  and  5.76  MHz,  which  reveal  distinct  BUM. 
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ERP  Dependence  of  BUM  -  Pump  =  5.67  MHz 


ERP  Dependence  of  DM  -  Pump  =  5.79  MHz 


Fig.  4.2.  Behavior  of  the  BUM  amplitude. 

Observations  made  by  the  Kodiak  radar.  Figure  4.3  a)  plots  time  series  of  the  SNR  of  reflected  as 
detected  by  the  Kodiak  radar.  We  chose  beam  9  directed  towards  HAARP’s  heated  region,  the 
frequency  of  the  probing  signal  is  about  12  MHz.  The  signal  reflected  by  the  artificial  ionospheric 
turbulence  is  centered  at  700  km  slant  range.  During  the  10*  cycle  at  03:12:00-03:14:30,  the  SNR 
was  strongly  reduced  since  Kodiak’s  beam  09  was  almost  out  of  the  HF  heating  spot.  This  was  a 
geometric  effect  of  the  angular  adjustments  of  the  beam  to  track  PRN  25  and  this  final  heating 
cycle  was  picked  up  by  beam  10.  The  HF  heating  was  switched  on  at  02:55  followed  by  the 
buildup  of  the  reflected  radar  signal  over  25  sec.  Then  a  strong  SNR  of  more  than  50  dB  was 
detected.The  reflection  increases  by  ~10  dB  around  02:57  then  it  dips  around  03:00  and  03:04  due 
to  the  decrease  in  the  ionospheric  plasma  density.  The  second  dip  was  reflected  in  the  ionogram  as 
discussed  above,  while  the  first  dip  occurred  during  the  time  not  covered  by  the  ionograms  at 
HAARP.  Between  03:1 1:30  and  03:13:30  the  SNR  decreases,  at  03:14:30  the  heater  was  switched 
off  and  the  reflection  gradually  disappeared. 

Figure  4.3  b)  shows  time  series  of  the  velocity  of  the  ionospheric  irregularities.  Starting  at 
02:55:00,  the  irregularities  quickly  accelerate  to  35-40  m/s,  peak  at  ~60  m/s  near  02:57:00- 
02:59:00,  then  gradually  reduce  to  ~30  m/s. 
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Fig.  4.3.  Kodiak  radar  observations 
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HAARP  signals  detected  at  UAS.  Monitoring  of  the  HAARP  HF  signals  was  carried  out  at  the 
Ukrainian  Antarctic  Station  (UAS)  “Academik  Vernadsky”.  Figures  4.4  a)  and  4.4  b)  show  the 
spectrogram  of  the  detected  signal  for  June  6  and  7,  2014,  while  figure  4.5  shows  the  intensity  of 
the  detected  signal  on  June  6*. 


June  06, 2014  (log  scale)  dB 


Time,  UT 

a) 


June  07,2014  (log  scale) 


Time.  UT 

b) 

Fig.  4.4.  Spectrograms  as  detected  in  Antarctica  at  the  “Akademik  Vernadsky”  station  for  June  6 

(a)  and  7  (b). 

The  8-th  (fh  =  5850  kHz)  and  9-th  (fh  =  5910  kHz)  heating  cycles,  indicated  by  black  arrows  at  the 
spectrograms,  were  contaminated  by  interference  signals  and  therefore  were  not  processed  in  Fig. 
4.6  which  reveals  the  signal-to-noise  ratio  of  the  signals  intensity.  Figure  4.5  shows  that  the 
intensity  of  detected  signal  strongly  depends  on  the  heating  frequency  fr.  In  fact,  for  fr  slightly 
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above  the  4^^  gyro  the  intensity  of  the  deteeted  signal  was  very  low,  hardly  above  the  noise  level. 
The  intensity  of  the  deteeted  signal  went  up  with  fh  and  peaks  at  5.79  MHz. 


Received  Signal  at  UAS  -  Verdanasky  -  06-Jun-2014  02:55:00  UTC 


480  600  720 

Elapsed  Time  (s) 


1200 


Fig.4.5.  Frequency  dependency  of  the  detected  signal  intensity. 


Fig.  4.6.  Signal  to  noise  ratio. 

Further,  Fig.  4.6  shows  how  the  signal  changes  with  the  intensity  of  the  HAARP  heater  at  the 
different  heating  frequencies  (the  blue  lines).  The  green  lines  correspond  to  4  times  increase  of  the 
heating  power,  the  linear  interpolation  is  shown  by  the  red  lines  in  both  Fig.  6.  Similar  results 
were  obtained  on  June  7  except  that  the  signal-to-noise  ratio  was  by  5-10  dB  less  than  that  on  June 
6  (see  Fig.  4.4  a  and  Fig.  4.4.  b). 
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The  spectral  width  as  a  function  of  the  heating  frequency  is  shown  in  Fig.4.7  for  each  of  100  s 
time  interval  which  correspond  to  the  heating  by  a  chosen  frequency. 


July  6.  2014  (Width  of  the  spectra) 


5670  5700  5730  5760  5790  5820  5850  5880  5910  5940 

F^,.  kHz 


Fig.  4.7.  Spectral  width  of  the  signals  detected  in  Antarctica. 

The  spectral  width  at  half  power  was  obtained  by  using  smoothed  spectral  curves.  It  was  found 
that  the  narrow  spectral  width  was  about  1  Hz  on  6/6/2014  and  increase  to  2  Hz  at  6/7/2014. 
Finally  we  obtain  the  S4  scintillation  index  over  100  s  time  interval  as  above.  It  is  presented  in  Fig. 
4.8.  as  a  function  of  the  heating  frequency  for  June  6  and  7  experiments. 


06.06.2014  02:55-03:15  UT 


0.6 


5670  5700  5730  5760  5790  5820  5850  5880  5910  5940 

Frequency,  kHz 


Fig.4.8.  S4  scintillation  index  of  intensity  of  the  detected  HAARP  signals. 

In  both  experiments  the  S4  index  peaked  at  fh=5.73  MHz  i.e.  when  the  SSS,  indicated  by  the 
BUM,  begin  to  be  suppressed  by  the  long  scale  striations  such  as  revealed  by  the  DM,  the  latter 
effectively  scatter  the  HF  waves. 
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Discussion .  Now  we  discuss  the  probing  of  the  AIT  by  HAARP  HF  waves  scattered  by  the  AIT 
into  the  ionospheric  waveguide.  The  radio  waves  were  detected  and  processed  at  Ukrainian 
Antarctic  Station  (UAS)  15.6  Mm  away  from  HAARP.  The  results  of  these  experiments  were 
presented  in  Section  3,  here  we  analyze  the  related  physical  processes.  Note  that  Zalizovski  et  al. 
[2009]  detected  HF  signals  from  EISCAT  facility  at  Ukrainian  Antarctic  Station  (UAS)  located 
16.3  Mm  from  the  heater.  The  HF  signal  revealed  two  different  components,  one  a  narrowband 
“mirror-reflected”  and  another  “scattered”  having  a  few  times  wider  spectrum.  It  was  mentioned 
by  Zalizovski  that  although  both  components  propagate  along  the  ionospheric  channel  the  mirror- 
reflected  signal  is  associated  with  the  side  lobe  radiation  of  the  transmitting  HF  antenna,  which 
bypasses  the  modified  volume,  while  the  scattered  signal  was  radiated  through  the  main  antenna 
beam  and  then  scattered  by  pump-induced  plasma  striations.  Notice  that  in  order  to  distinguish 
which  of  the  two  mechanisms  dominates  one  needs  to  sweep  the  heating  frequency  around  the 
multiple  gyro  resonance.  While  the  sidelobe  radiation  is  not  sensitive  to  the  frequency  variations, 
the  spectrum  of  the  pumped  striations,  so,  the  scattered  signal  strongly  depends  on  the  fr  choice. 

Figure  4.9  shows  amplitudes  of  PSDs  obtained  from  the  SEE  spectrum  combined  with  the 
intensity  of  the  HAARP  signal  measured  at  UAS  taken  from  the  Doppler  spectrum  (Fig.4.4) 
versus  the  heating  frequency. 

Normalized  SEE  and  Recieved  Power  at  UAS 


Fig.4.9.  Amplitudes  of  PSDs  of  SEE  and  intensity  of  the  HAARP  signal  measured  at 

UAS. 

The  red  and  blue  lines  in  Fig  4.9  show  the  amplitude  of  the  BUM  and  DM  respectively, 
both  normalized  by  their  peak  values.  The  black  line  shows  the  intensity  of  the  detected  HAARP 
signal  first  smoothed  by  moving  average  (red  trace  in  Fig  4.5)  and  then  normalized  by  its  peak 
value. 
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Note  that  the  Doppler  broadening  deteeted  on  June  6  reveals  that  AfD=1.15-1.40  Hz.  This 

cAf 

corresponds  to  the  velocity  of  ionospheric  irregularities  v  = - ^  =  30 -35  m/s,  and  it 


2/. 


IS 


consistent  with  the  velocity  found  by  Kodiak  in  figure  4.3.  It  was  found  that  the  intensity  of  the 


signal  detected  at  UAS  during  the  experiments  on  June  6  and  7  strongly  depends  on  the  heating 
frequency  fh.  Figure  9  reveals  that  on  06/06  when  the  BUM  is  dominant  in  the  SEE  spectrum  it 
suppresses  the  HF  signal.  The  stronger  is  BUM  the  weaker  is  the  HF  signal.  To  the  contrary,  when 
the  DM  dominates  in  the  SEE  spectrum  intensity  of  the  detected  signal  went  up  with  fh  and  peaks 
at  5.79  MHz.  Thus  the  DM  and  HF  peaks  occur  at  the  same  frequency. 

Furthermore,  it  is  known  that  the  BUM  is  associated  with  the  pumping  of  10  cm  transverse 
scale  striations,  often  called  super  small  striations  (SSS),  while  the  DM  is  associated  with  the  7-30 
m  scale  striations  [Norin  et  ah,  2008].  Since  SSS  inefficiently  scatter  the  HF  waves  compared  to 
the  decameter  striations  related  to  the  DM,  this  implies  that  the  HAARP  signal  detected  at  UAS 
could  be  due  to  the  scattering  of  the  HF  radiation  having  the  half  wavelength  of  25  m  by  the 
decameter  size  artificial  striations  into  the  ionospheric  channel. 

The  analysis  shows  that  origin  of  the  detected  signal  at  the  “Akademik  Vernadsky”  station 
is  scattering  of  the  HAARP ’s  HF  radiation  off  of  artificially  pumped  striations  and  into  the 
ionospheric  waveguide.  In  this  process,  mirror  reflection  does  not  appear  to  play  an  important 
role,  as  indicated  by  the  low  amplitude  of  received  signal  at  5.67  and  5.70  MHz.  Since  the 
amplitude  of  the  signal  varied  nonlinearly  with  the  HF  power  at  5.76,  5.79  and  5.82  MHz,  the 
signal  cannot  be  induced  by  sidelobe  radiation  which  is  known  to  be  linear  in  the  radiated  power. 
Still  the  question  remains  why  the  observed  spectrum  is  relatively  narrow  about  1-2  Hz  width, 
which  is  similar  to  that  observed  by  Zelizovski  et  al.  [2009]  of  the  mirror-reflected  signal 
component  while  the  scattered  signal  was  much  broader  up  to  10-20  Hz. 
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5.  Conclusions 

The  Project  was  aimed  at  experimentally  investigating  the  possibility  of  exciting  the  ionospheric 
interlayer  duct  channel  using  powerful  radiations  from  the  heaters  HAARP  (Alaska,  USA)  and 
EISCAT  (Tromso,  Norway),  as  well  as  from  HF  broadcasting  stations  RWM  (Moscow,  Russia) 
and  CHU  (Ottawa,  Canada).  Major  attention  was  paid  to  analyzing  the  possibility  of  exciting  the 
interlayer  ionospheric  waveguide  which  support  super-long  range  HF  propagation  with  a  small 
amount  of  attenuation.  To  monitor  the  radiation  a  compact-size  receiving  complex  was  developed 
which  is  capable  of  measuring  the  signal  intensity  and  spectral  characteristics  in  an  off-line 
automatic  mode.  Two  facilities  of  the  kind  have  been  constructed  in  the  course  of  the  Project 
fulfillment.  One  of  these  was  deployed  in  Ukraine  at  the  Low-frequency  Observatory  of  the  IRA 
NASU  (Martova  village,  Kharkov  region)  in  2012,  while  another  was  installed  at  the  Ukrainian 
Antarctic  station  “Akademik  Vernadsky”  in  the  Antarctic  in  2013.  In  addition,  similar  receiving 
complexes  were  involved  during  the  Project-related  measuring  campaigns,  which  had  been 
constructed  by  the  Project  participants  under  other  research  projects  and  put  into  operation  in  the 
Northern  Scandinavia  (Tromso,  Norway),  on  the  Svalbard  Island  (Spitsbergen  Archipelago)  and 
Equatorial  Africa  (Lagos,  Nigeria).  Such  a  chain  of  the  receiving  and  transmitting  sites  made  it 
possible  to  simultaneously  investigate  probe  signal  propagation  along  radio  paths  of  various 
lengths,  specifically,  quasi-vertical  (less  than  100  km),  middle  (up  to  3000  km)  and  super-long 
range  (from  10000  to  16000  km).  In  all,  about  100  hours  have  been  spent  for  observations  of  the 
radiation  from  powerful  heaters  (mostly,  of  the  EISCAT  heater)  and  more  than  3000  hours  for 
monitoring  signals  from  broadcast  radios.  In  a  number  of  cases  the  signal  strengthening  was 
detected  for  the  super-long  range  radio  links  (Alaska-Antarctica  and  Northern  Scandinavia- 
Antarctica)  which  effect  can  be  regarded  as  a  result  of  the  waveguide  propagation.  To  explain  the 
mechanisms  of  the  waveguide  excitation  and  extraction  of  the  radio  wave  energy  from  it,  the 
effect  of  aspect-sensitive  scattering  of  signals  by  natural  and  artificial  field-aligned  irregularities 
of  the  ionospheric  plasma  was  considered  as  the  basic  one.  A  pioneering  feature  of  the  developed 
theoretical  model  is  the  account  of  the  regular  ionospheric  refraction.  The  aspect-sensitive 
contours  in  the  ionosphere  and  on  the  Earth’s  surface  have  been  calculated  for  all  the  transmitting 
and  receiving  sites  for  the  current  ionospheric  conditions.  During  the  measurement  campaigns 
with  the  use  of  radiations  from  the  powerful  heaters  two  other  effects  have  been  occasionally 
revealed,  which  are  the  combination  (Brillouin)  signal  scattering  by  the  artificially  stimulated 
plasma  turbulence  and  excitation  of  the  second  harmonic  of  the  powerful  radiation  in  the  HF- 
modified  ionosphere.  As  planned,  a  software  package  has  been  developed  for  the  remote  control 
of  the  receiving  complexes  and  visual  representation  of  the  measurement  results  in  real-time 
through  the  internet.  It  should  be  noted  that  throughout  2013  the  HAARP  heater  was  out  of 
operation  because  of  the  lack  of  the  operational  cost.  In  view  of  this,  all  the  measurement 
campaigns  on  monitoring  the  powerful  transmitter  signals  in  the  Antarctic  have  been  carried  out 
using  the  EISCAT  heater  under  conditions  which  were  not  always  optimum  for  the  ionospheric 
waveguide  excitation.  For  the  same  reason  the  receiving  complex  for  the  Antarctic  has  been 
provisionally  deployed  at  the  Ukrainian  Antarctic  station  “Akademik  Vernadsky”.  The  most 
fruitful  experiments  were  performed  in  a  frame  of  BRIOCHE  heating  campaign  in  June  2014.  The 
detail  description  of  this  campaign  is  included  in  the  Final  Report  (Chapter  4).  All  results  the 
previous  EISCAT  heating  campaigns  were  included  into  the  Interim  Project  report. 

The  results  obtained  during  the  Project  fulfillment  have  been  published  in  fours  scientific 
papers  and  reported  at  several  international  meetings  in  the  USA,  Puerto-Rico  and  Ukraine. 
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